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We study how decoherence rules the quantum-classical transition of the Kicked Harmonic Oscil- 
lator (KHO). When the amplitude of the kick is changed the system presents a classical dynamics 
that range from regular to a strong chaotic behavior. We show that for regular and mixed classical 
dynamics, and in the presence of noise, the distance between the classical and the quantum phase 
space distributions is proportional to a single parameter x = Kti^g /4D'^''^ which relates the effective 
Planck constant hcs, the kick amplitude K and the diffusion constant D. This is valid when x < !> 
a case that is always attainable in the semiclassical regime independently of the value of the strength 
of noise given by D. Our results extend a recent study performed in the chaotic regime. 
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I. INTRODUCTION 



The complete description of the emergence of the clas- 
sical world within the quantum theory is a fundamental 
problem that have attracted a lot of attention since the 
beginning of quantum mechanics. Although many ad- 
vances have been made, there are still open questions. An 
important step forward was done with the understanding 
that the formalism of phase space distribution functions 
is the appropriate framework to study the quantum to 
classical transition Among all quantum phase-space 
distributions, the Wigner function |3| has the advantage 
that identifies the nonlinearities of the system as the re- 
sponsible for the separation between the quantum and 
the classical evolution. Indeed, while for linear systems 
the Wigner distribution obey the Liouville equation of 
the classical distribution, the nonlinearities add terms to 
this equation that eventually set the two distributions 
apart Q, even for initial states that are classically al- 
lowed. 

The time for which classical and quantum distributions 
starts to differ, called breaking or separation time, can be 
very short for classically chaotic systems. This is because 
the nonlinearities of the system are reached very quickly 
due to the exponential stretching of the distributions. 
Hence, the breaking time scale as ln(l/?icff), where the 
effective Planck constant hcs = h/S is the semiclassical 
jarameter (with S a typical action of the system) 0, |^ 



parameter (w: 



An important and subtle problem is that the break- 
ing time for macroscopic systems (hcs <C 1) can still be 
small compared to a typical evolution time due to its 
logarithmic behavior, allowing quantum effects that are 
not observed in the classical world llj. This paradoxi- 
cal situation is explained by the coupling of the system 
with the environment, which leads to the elimination of 



the quantum signatures, providing the reconciliation of 
theory and observation 0,^]. At this respect, an open 
question is how the effects of the environment affect the 
logarithmic law of the breaking time in systems with a 
classical chaotic dynamics. More precisely, we can say 
that it is still unknown the exact relation between the 
strength of noise and the effective Planck constant in or- 
der to obtain the correct classical limit ■ 

This problem has started to be answer in a recent 
work T3| where the authors study the effect of a purely 
diffusive environment over an specific model given by 
the Kicked Harmonic Oscillator (KHO). They show that 
the differences between the classical and quantum phase 
evolution is proportional to a single parameter x = 
Khlg/4D^^^ , that combine HeS, the diffusion coefficient 
D and the kicking strength K, who control the macro- 
scopicity, the noise and the chaotic behavior respectively. 
These quantum-classical differences were estimated with 
Vn {n is the number of kicks) which is the integral over 
the whole phase space of the modulus of the difference be- 
tween the quantum and the classical distributions. Thus, 
it was shown that, when the classical dynamics is chaotic, 
in the semiclassical limit hcs ^ 1 and when x < 1, the 
single parameter x controls the quantum-classical tran- 
sition of the KHO. This result, confirm the conjecture 
presented in Ref. that in the presence of noise a sin- 
gle parameter controls the quantum-classical transition 
of classically chaotic systems. 

In this paper we analyze the case when the classical dy- 
namics of the KHO is regular or mixed. In the absence 
of a reservoir, we show that the separation between the 
quantum and the classical evolution can be characterized 
through the behavior of I?„ as a function of time. When 
the initial distribution explore regions of regular classical 
dynamics, the behavior of allows the identification 
of a potential scale-law f^or the separation time. 



However, this scale-law is not valid when the degree of 
mixing increase. In the presence of a purely diffusive 
reservoir, we show that the quantum-classical differences 
measured by still scale with the parameter x in the 
semiclassical limit and when x < 1. Hence, the con- 
cept of separation time is not meaningful anymore, since 
the quantum an the classical distributions remain close 
throughout the evolution. Our result suggest that the 
conjecture in can be extended independently of the 
underlying classical dynamics. 

The paper is organized as follows. Sectionnjis devoted 
to review the quantum and the classical versions of the 
KHO and also to explain the formalism for the evolution 
of the Wigner function and the classical distribution as a 
map in phase space. In Section lTTll we introduce our mea- 
sure of the quantum-classical differences I?„, and study 
its behavior as a function of time without the presence 
of a reservoir. We investigate the scale- law with the ef- 
fective Planck constant of the separation time when the 
initial distribution explore regular or mixed regions of the 
phase space. The effects of decoherence over when 
the KHO is in contact with a purely diffusive reservoir, 
is studied in Section llVl We end with some final remarks 
m Section El 

II. MODEL: THE KICKED HARMONIC 
OSCILLATOR 

A. The quantum version 

The model used in our calculations is the KHO, a par- 
ticle of mass rn in a harmonic potential subjected to a pe- 
riodically applied position-dependent delta pulses. The 
quantum Hamiltonian is defined as: 

H = + -mi/^Q^ + A cos(kQ) V S(t - nr), (1) 

2m 2 

n=0 

where i/ is the oscillator frequency, A the amplitude of the 
kicks and r the interval between two consecutive kicks. 
k sets the periodicity of the position dependent kicking 
potential. The physical realization of this Hamiltonian 
is within the context of trapped ions. It was shown in 
Ref. 17], that it describes the center-of-mass dynamics 
of an ion in a one-dimensional trap submitted to a se- 
quence of standing-wave laser-pulses, off-resonance with 
a transition between the electronic ground state and an- 
other internal state. The wave number k in Eq. ^ is the 
projection along the trap axis of the corresponding wave 
vectors with equal modulus | fc | of two opposite propagat- 
ing pulses with oblique incidence measured by an angle 
0, i.e. k = 2|fc|cos(6'). 

It is possible to define an effective Planck constant 
hcs if we work with dimensionless quantities q = kQ, 
p — kP/mv, and K = k^A/mv, so that [q,p] = 2iif' = 
ihcs, with 7] = fcAQo = k^h/2mv., and AQq being the 
width of the ground state of the harmonic oscillator. The 




FIG. 1: Stroboscopic phase space of the KHO for different 
values of the kicked amplitude K. In (a) K = 0.5, (b) K = 
1.0, (c) K — 1.5 and in (d) K — 2.0. The phase space points 
(gOiPo) = (0,0), (0,1.1) and (0.7,0), marked with symbols 
(•), are the centers of the initial coherent states used in our 
numerical simulations. 



dimensionless parameter 77 is the so-called Lamb-Dicke 
parameter [l8l |. which measures the ratio between the 
ground state width and the wavelength A = 2n/k that 
sets the scale of the non-linearity of the Hamiltonian. In 
experiments with trapped ions, the classical limit ry — > 
(and thus = 2r]^ — > 0) can be approximated sim- 
ply by changing the angle 9 of incidence of the incoming 
pulses, or by increasing the trap frequency v. 



B. The classical dynamics. 

The classical Hamiltonian is obtained from Eq. ^ sim- 
ply by replacing the quantum operators Q and P by the 
phase space variables Q and P respectively. In an anal- 
ogous way we can define the dimensionless phase space 
coordinates x = {q,p) = {kQ, kP/mv). Then, the time- 
dependent classical evolution can be described as a com- 
position of a discrete map corresponding to the kick, plus 
a rotation in phase space: x„_|_i = R o K(x„) where 
x„ = (g„,p„) are the coordinates before kick n (where 
the first kick corresponds to n = 0). The kick operation 
K is defined as 

qt =1n, Pt =Pn + K sm(qn) , (2) 



3 



and the phase space rotation R is given by: 

<?«+! = cos(i/t) q+ + sin(i/r)p+ , 
Pn+i = - sin(i^r) q+ + cos(iyT) p+ 



(3) 



In close relation with Ref. [T^ , we also consider the case 
where the interval between kicks r and the period of the 
harmonic oscillator T = I-k jv are related by r = T/6. 

The phase space of the system is unbounded and cor- 
responds to a mixed dynamics exhibiting stable islands 
surrounded by a stochastic web along which the particle 
diffuses. The thickness of the web is controlled by the 
chaoticity parameter K = k'^A/mv which is the dimen- 
sionless kicked amplitude. Hence, these chaotic regions 
broadens (shrinks) as K increases (decreases). In the 
case of r = T/6 considered, the stochastic web displays 
crystal hexagonal symmetry In Fig^we show the 

stroboscopic phase space around the origin for different 
values of K. For K < Kc ~ 1.1547 the origin is an 
elliptic fixed point while for K > Kc the system under- 
goes a bifurcation where the origin become an hyperbolic 
fixed point and two neighbor elliptic fixed points arise. 
The separation between these two elliptic fixed points in- 
creases with K. As we can see in Fig^a), for K = 0.5 
the dynamics of the system around the origin is essen- 
tially regular with a series of concentric tori around the 
elliptic fixed point. For K — 1 [Fig^Jb)], some tori were 
destroyed, and tiny chaotic regions appear (invisible to 
the eye). Mixed dynamics star to manifest most appre- 
ciable for K = 1.5 while ior K = 2 most of the area 
around the origin shows chaotic behavior [see Fig ^^c) 
and (d)]. 



C. The evolution of phase space distributions 
without reservoir. 

We analyze the quantum-classical transition in phase 
space comparing the evolution of the Wigner function 2] 
and the corresponding classical distribution. The Wigner 
function, defined as 



W(x) = dy{q + y/2 



P 



Airrf 



q - 2//2 ) cxp <^ -i-T-^ 



py 

'V, 
(4) 

with p the density operator of the quantum state, 
have a highlight position within all the quantum quasi- 
probability distribution in phase space because it is the 
only one to yield the correct marginal probabilities, thus 
satisfying an important property of classical distributions 
|l9j| . Furthermore, quantum effects are more pronounced, 
as compared to other bounded distributions. Finally, new 
methods have been developed for the direct measurement 
of the Wigner function, which are adequate for the de- 
scription of evolving systems . 

Without the action of a reservoir the unitary evolution 
of the Wigner function of the KHO can be written as 



iy„+i(x) = J dx' L(x^,x') P^„(x') 



(5) 



where Wn+i and Wn are the Wigner functions imme- 
diately before the kicks n + 1 and n respectively. The 
integration is over all the unbounded phase space, and 
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(6) 



is the quantum propagator of one kick plus the harmonic 
evolution between the consecutive kicks which is given 

by 



x«^[g^'(x),p^(x)]^R-i(x) 



(7) 



i.e. the phase space coordinates rotated with the inverse 
transformation of Eq. ^ . The Liouville evolution of the 
classical distribution 14^^' (x) can also be written in the 
form of Eq. (jSJ with the classical propagator 

L^' (x^ , x') ^S[p'- + K sin(g')] 5 (g^ - q') ■ (8) 

In the classical limit 77 ^ the classical propagator is 
formally recovered from the quantum one 5] . Indeed, in 
the semiclassical regime (77 ^ 1) the integral in Eq.© 
can be estimated employing stationary-phase techniques 
yielding. 



L(x«,x') 



1 



Ai 



|6|i/3 



sign(6) 
' |6|V3 



Ksm{q')) j X 
(9) 



where b = {rj'^K/2) sm{q') and Ai(a;) is the Airy function. 

When 77 — > we can use that Ai{— y/e}/e S{y) to get 
the classical Liouville propagator in Eq.®. However, no 
matter how small we can do the Lamb-Dicke parameter rj 
in a physical realization of the KHO, this will always take 
a finite value and then the presence of an Airy function in 
the quantum propagator will eventually set the Wigner 
function and the classical distribution apart. 



III. THE QUANTUM-CLASSICAL 
SEPARATION WITHOUT RESERVOIR: 
SCALE-LAW FOR THE SEPARATION TIME. 

The Wigner distribution of a linear system (quadratic 
Hamiltonian) satisfy the Liouville equation for the clas- 
sical phase-space probability distribution. Thus, if we 
start with a Wigner function of a classically allowed prob- 
ability distribution (the only possibility is a minimum 
uncertainty initial Gaussian wave packet |22j'). the lin- 
ear dynamics evolution will preserve its Gaussian shape. 
When the system have non-linearities, the initial wave 
packet will deform and eventually delocalized during the 
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evolution. As a consequence, quantum interference ef- 
fects appear between different pieces of the wave packet, 
which remains coherent throughout the unitary evolu- 
tion. These quantum interferences effects stem from the 
extra terms that the non-linearities of the system add 
to the Liouville equation in order to obtain the correct 
equation of motion of the Wigner function 0,[23- Equiv- 
alently, we can say that these interferences are originated 
by the non-local nature of the Wigner propagator (see for 
example Eq.©] in comparison with the local nature of 
the Liouville propagator of Eq.©, that persist even in 
the semiclassical regime [see Eq.®]- Hence, they are 
the responsible for the separation of the Wigner function 
from the classical distribution and thus of the quantum- 
classical correspondence breakdown. 

The quantum-classical differences start to become im- 
portant when the initial distribution spreads over dis- 
tances of the order of the characteristic scale of the 
nonlinearities of the system, and it is well known this 
argumentation in order to estimate the time scale for 
the separation in the case of classically chaotic systems 
[iLI2»24j- Indeed, as the initial Gaussian Wigner func- 
tion is well localized, it will initially evolve following 
classical trajectories in phase space, therefore the wave 
packet will stretch exponentially fast along the unstable 
manifold and will shrinks also exponentially fast along 
the stable one in order to conserve the phase space vol- 
ume. This is also accompanied by the characteristic fold- 
ing of classical trajectories that forces the wave packet to 
develops quantum interference between different pieces 
which are coherently related. As a result of this behav- 
ior, it is found a logarithmic scale law with the effective 
Planck constant hcs of the separation time in systems 
with classical chaotic dynamics. For the KHO, the esti- 
mation of the separation time tg gives 

(10) 

where A is the local expansion coefficient and fioff = 277^ . 
A detailed estimation of ts can be found in Ref. 0, 0| 
and in , where it is also shown a numerical test that in- 
dicates that this is also the time when quantum and clas- 
sical expectation values start to differ from each other. 
The expansion coefficient A corresponds to the Lyapunov 
exponent when an average over initial conditions in the 
chaotic regions is considered 10, 15] . 

In the other extreme situation, when the classical dy- 
namics of the system is essentially regular, the scaling 
law of the quantum-classical separation time with hos 
was less studied. The logarithmic dependence of tg in 
the chaotic case is due to the exponentially divergence of 
classical trajectories, but in the integrable case the diver- 
gence is typically polynomial 0,l23- Therefore, follow- 
ing a similar type of argumentation that in the chaotic 
cas e, the separation time in regular systems would follow 
a potential scale-law with h^g. In the case of the KHO, 
if the initial distribution explores regions of essentially 
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FIG. 2: Separation D„ between quantum and classical distri- 
bution as a function of the number of kicks for K — 0.5. The 
full lines joint all the values of Dn calculated for the same 
value of the Lamb-Dicke parameter rj. In (a) the initial dis- 
tribution is a coherent state centered at the phase space point 
(gO,pO) — (0, 1.1) and the curves are for the decreasing values 
of 7? = 0.5, 0.25, 0.125, 0.0625 and 0.03125. The result for the 
initial condition (gO,pO) — (0.7,0) is equivalent (not plotted 
here). In b) the initial coherent state is centered at the origin 
and the curves are for 77 = 0.5,0.25,0.125 and 0.0625. The 
number of kicks n are rescaled by a function n*{rj) = ri'" 
where we optimize the values of a in order to obtain the best 
collapsing of all the curves. The insets show the curves with- 
out the rescahng over the time n. 

regular dynamics, we would have 

= UsT « r/r?" , (11) 

where numerical constants are forbidden. Indications of 
the validity of Eq. l(TT|l for if <C 1 are given in Ref. . 
In what follows we show a way to estimate the amount 
of quantum effects that lead to the separation between 
the Wigner function and the classical distribution, and 
how the scaling law for the separation time in Ec ljlll) 
manifests in our formalism. 

In Ref. jl5l |. it was introduced the time dependent 
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FIG. 3: Idem as FigEl but for K = 1.0. In a) the ini- 
tial coherent state centered at (gO,pO) — (0,1.1) and in b) 
(gO,pO) — (0.7,0). The different curves in each plot corre- 
sponds to 77 = 0.5, 0.25, 0.125, 0.0625 and 0.03125. The num- 
ber of kicks n are rescaled by a function n*{rj) — r;~" where 
we optimize the values of a in order to obtain the best col- 
lapsing of all the curves. 

quantity (t = nr), 

I?„ = J dx |M^„(x) - W,fi^)\ , (12) 

to quantify the quantum effects that lead to the separa- 
tion between the Wigner function and the classical distri- 
bution. Here, Wn and W^' are respectively the Wigner 
function and the classical distribution immediately be- 
fore the kick n, both normalized to the unity. Starting 
with the same initial distribution it is expected that dur- 
ing the unitary evolution the appearance of quantum- 
classical differences make I?„ to grow and eventually sat- 
urate around some finite value. This behavior is shown 
in Fig|21 and 13 where we have plotted I?„ as a function of 
the scaled number of kicks (see below) for different val- 
ues of the semiclassical parameter 77 and the amplitude of 
the kick K. The initial distribution is always a coherent 
state with width Aq{0) = Ap{Q) = rj. 



In Fig|21the amplitude of the kick is -ft' = 0.5, so we 
study the case with regular classical dynamics as the dif- 
ferent initial coherent states considered, during the in- 
terval of times monitored, spread over the regular region 
around the origin (see Fig. ^ 2]and[5ll. We can see in 
Fig. 121 (a) that when we start with an initial wave packet 
outside the origin, the quantity I?„ grows up to values 
larger that one, for all values of the semiclassical param- 
eter T] [2^ . This clearly indicates the separation between 
the Wigner and the classical distributions because a mea- 
sure of their differences reach values of the order of their 
normalization. However, it is important to note that no 
matter which value of the quantum differences, measured 
by Vn, we choose in the plot of Fig.|21(a) to define the sep- 
aration between the Wigner and classical distributions: 
it will be reached at some early time of the evolution in- 
dependently of how small is the semiclassical parameter 
77. Thus, the separation time defined in this way should 
scale as in Ea. (|ll|l . This is showed in the plots by the 
rescaling of the number of kicks by a function of the form 
n*{r]) = ri~°' that allows the approximately collapse of all 
the curves. In Fig.|2|it is shown the case when the kick- 
ing amplitude is K = I. We can see a similar behavior 
although the collapsing of all the curves do not work so 
well. This is an indication of the departure of the separa- 
tion time from the scale-law in Ea. (|ll|l when the degree 
of mixing in the classical dynamics grows. It is important 
to note that such scaling was not observed for the case 
when the kicking amplitude is if = 1.5. 

In Fig. 12 (b) we show the case when the initial coher- 
ent state is centered at the origin of phase space when 
K — 0.5. Vn grows up to values around 0.25 and then 
saturate for all semiclassical parameters 77. In this case, 
Vn grow less than when the initial coherent state is out- 
side the origin because now the quantum and classical 
distributions remains well localized for very long times 
and thus the quantum-classical differences are less pro- 
nounced (Fig. 01. Nevertheless, this differences do not 
decrease for small values of 77, showing that the sepa- 
ration between the classical and the quantum evolution 
occurs in the semiclassical regime. The fluctuations ob- 
served in Fig. 121 (a) and (b) have a period of 77 « 6 kicks 
i.e. approximately the period of the harmonic evolution 
T = 6r. This is the time that most parts of the quantum 
and classical distribution take to complete a spin along 
the elliptical structure of phase space around the origin 
for K = 0.5 (see Fig.lU. 

It is interested to compare in Fig. 21 and Fig. [SI the 
quantum-classical differences that leads to the separation 
between the quantum and the classical evolution when 
the amplitude of the kick is K — 0.5. Indeed, when 
the initial coherent state is centered outside the origin, 
the quantum and classical distributions spread out and 
as a consequence the Wigner function develops a quan- 
tum interference pattern throughout the evolution that 
clearly indicates its separation from the classical distri- 
bution (see Fig. 2} . Now, when the initial coherent state 
is centered at the origin of phase space, which is an el- 
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FIG. 4: Density plots of the Wigner (right) and classical dis- 
tributions (left) that result from the evolution with the KHO 
Hamiltonian with the kick amplitude K = 0.5 from an initial 
coherent state centered in the phase space point (go,Po) ~ 
(0.0, 1.1) and with a width Aq{0) = Ap{0) = 77 = 0.25. In (a) 
and (b) immediately before the kick n = 8 and in (c) and (d) 
immediately before the kick n = 18. 



liptical fixed point for K = 0.5, the quantum and the 
classical distributions remain well localized for very long 
times (Fig. (Sjl. In this case, the classical distribution 
develops the typical classical structure of "whorls" asso- 
ciated with an elliptic fixed point 0, |^ . It is evident 
from the figure that the quantum evolution smooths out 
this details of the classical distribution and this is the 
origin of the quantum-classical differences that leads to 
the separation. 



IV. DECOHERENCE: THE EFFECT OF A 
DIFFUSIVE ENVIRONMENT 



Here, we study how the decoherence affects the 
quantum-classical differences given by the quantity I?„, 
when the initial distribution explore phase space regions 
of regular or mixed classical dynamics. Our approach is 
in line with the one developed in Ref. 15], where the 
case of chaotic dynamics was studied. Therefore, we also 
coupled the KHO to a thermal reservoir with average 
population n, in the Markovian and weak coupling limit, 
and we consider the purely diffusive regime '15], i.e. the 
high temperature regime n — > 00 together with F ^ 
but when nT is constant (F is the dissipation rate). In 
this limit, the action of the reservoir over the Wigner 



FIG. 5: Density plots of the Wigner (right) and classical dis- 
tributions (left) that result from the evolution with the KHO 
Hamiltonian with the kick amplitude K = 0.5 from an initial 
coherent state centered at the origin of phase space and with 
a width Aq{0) = Ap(0) = rj ^ 0.25. In (a) and (b) immedi- 
ately before the kick n — 60 and in (c) and (d) immediately 
before the kick n — 300. 



function is described by the Fokker-Planck equation 



dW 



= F 



reservoir 



(13) 



with F = nVif' . For the complete evolution we have 
to add to right-hand side (RHS) of Ea. (|13|l the unitary 
evolution given by the Hamiltonian in Eq. ^ . The action 
of a pure diffusive reservoir over a classical distribution 
is also described by Ea. H13|) if we identified F with the 
classical diffusion constant. In this case, we must add to 
the RHS of Eg. lfT^ the corresponding Liouville term of 
the classical dynamics to obtain the complete evolution. 

Our findings about the behavior of the quantity 2?„ in 
the presence of decoherence are a consequence of the ef- 
fects of the reservoir over the one kick propagator of the 
Wigner function, that bring it close to the classical prop- 
agator of one kick. Hence, we first summarize the steps 
in Ref. that allows to write the one kick propagator 
of the Wigner function as a function of the one kick prop- 
agator of the classical distribution under the action of a 
purely diffusive reservoir and in the semiclassical regime 
(?7< 1). 

The solution of the differential equations for the evo- 
lutions of both, the classical and the quantum distribu- 
tions, in the interval of time between two consecutive 
kicks, and in the presence of the purely diffusive reser- 
voir, can be written as in Eq.Q where we have to replace 
the one kick propagator L(x,x') of Eq.® or ||SJ| by the 
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FIG. 6: Renormalized distances between quantum and classical distributions, with diffusion. All the curves were generated with 
the initial coherent state centered in the phase-space point (go,Po) ~ (0, 1.1) (see Fig. 0) with the kicked amplitude K = 0.5 
in (a) and K = 1.5 in (b). The eleven different curves correspond to x = Kri*/D^^^ = 6.2 x 10"^ 4.3 x 10"^, 1.5 x 10"^, 1.1 x 
10-^ 7.6 X 10-^ 3.9 X 10-^ 1.4 X 10"^, 6.8 x lO"*, 4.8 x 10-"*, 2.4 x 10"'', 4.3 x IQ-'^ in (a) and x = 6.5 x 10-^ 4.5 x 10-^ 3.3 x 
10~^2.3x 10"^ 1.2 X 10"^4.1 x 10"^2.1 X 10~^ l.lx 10~^7.2 x 10"", 1.3 x 10"", 4.5 X 10"^ in (b), where we used the values 
of diffusion constant D = 0.1, 0.05, 0.01, 0.005, 0.001, 0.0005 and the Lamb-Dicke parameter rj = 0.25, 0.125, 0.0625, 0.03125. 
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FIG. 7: The maximum distance between quantum and classical distributions given by the value of the first peak of T>n at 
n = ripeak as a function of x. The initial distribution is centered at xo = (0, 1.1) in (a) and xo = (0.7, 0) in (b). The marks (o) are 
for K = 0.5, the (■) for K = 1.0 and ( A) for K = 1.5. The region of linear behavior of I'npeak ^ function of x is highlighted by 
the dashed line with unit value of the slope. The inset is the result borrow from 15j when the initial distribution explore a chaotic 
region in phase space of the classic KHO. We used the value of the Lamb-Dicke parameters rj — 0.5, 0.25, 0.125, 0.0625, 0.03125 
and the diffusion constant D = 5 x 10"^ 1 x 10"^ 5 x 10"^ 1 x 10"'', 5 x 10"'', 1 x 10"^ 5 x 10"^ 1 x 10~^ 5 x 10~^ 10"^ 



smoothed one 

here D = Fr is the diffusion constant between two con- 
secutive kicks, and x-'^ is the result of the harmonic evolu- 
tion given in Eq.©. Then, the classical smoothed prop- 
agator for the classical distribution function is 

i^'(xV)-^^, (15) 



where y — {p' ^ + K sin g') j^^pD and x = 
(g' — g^) /2\fD. In the same way, we performed the in- 
tegration in Eq. (|14|l over the the propagator in Eq. (jHJ to 
obtain the smoothed quantum propagator: 



L(x^x')^ J^e-'^e^[''™'(''^^'"(^^-''(^"-^'^U(g«-g') ■ (16) 



8 



The presence of the Gaussian factor exp(— Z)/x^/?7**) 
in the integrand allows to obtain an approximation in 
terms of the classical smoothed propagator given in 
Ea. (|15|l . Indeed, when the width of this Gaussian is 
small, if'l^lD <^ 1, the /i's that effectively contribute 
to the integration are those close to the origin. So, it can 
be used sin(/i) w /i — in the phase of the integrand. 
Moreover, if x = Krj'^/ D^^^ <^ 1, the term with fi^ in the 



ie 



phase is small, so it also can be used e*'-^"'"'^-' w e'" + iS e 
and then the /i-integration can be performed. Thus, it 
is obtained the following approximation to the smoothed 
quantum propagator: 



L(x«, x') « i^'(x«, x') [1 + X sin(g')/(2/)] 



(17) 



where / {y) = 1/4 {y - 2y^/3) . Since | / iy)exp {-y^) \ < 
0.081, Eq. H17|l is valid under the less restrictive condi- 
tion X ^ 1- It is important to note that the smoothed 
quantum propagator (Eq.l^T))) is obtained without any 
assumption of the underlying classical dynamics of the 
system. The only conditions that have to be fulfilled are: 
i) £> > 77'' and ii) D > (Kri'^)^^^. These two conditions 
can be casted in D > {Krj'^)^/^ ^ •q'^, which is always 
attainable in the semiclassical regime (77 ^ 1), indepen- 
dently of the value of D. 

The approximation of the quantum propagator 
[Ea. H17|) ] allows us to write the Wigner function immedi- 
ately before the kick n as. 



W„(x) = W^\^) + X V G,(x) + 0{x^) 



(18) 



where the tilde indicates the action of the purely diffusive 
reservoir and we define the phase space function, 

G, (x) = j dx„L^'(x«,x„) 

X j dxj sin((7j)/(%) Z"'(xf ,xj_i) X 



dxo Z='(xf,xo) Wo(xo) 



(19) 



Ea. (|18|l is obtained iterating the evolution given by 
Eq. (|SJ with the smoothed Wigner propagator in Eq. H17() 
from the initial condition M/o(xo)i a-rid keeping only terms 
of first order in x- Hence, when x 1 '^s get for the sep- 
aration between the classical and quantum distributions. 



dx 



(20) 



The scaling law proportional to x is confirmed in Fig. El 
where we displayed I'n/x ^ function of the number of 



kicks n for K — 0.5 and K = 1.5 and a wide range of 
values of rj and D. We can see that for order of magni- 
tude different values of x, all the curves I?„/x-vs.n fit in 
the same scale. So, independently of the value of D, in 
the semiclassical limit we have D » {Ki^'^ f/^ » t?^, and 
thus the separation between the quantum and the classi- 
cal distributions goes down with 77^, becoming arbitrarily 
small. The curves in Fig. were generated for an initial 
coherent state whose Gaussian Wigner function, W^o(xo), 
is centered in the phase space point Xq = (0, 1.1) (Fig.0 
but we get similar results when Wo(xo) is centered at 
Xq = (0.7, 0). It is worth to remember that the initial dis- 
tributions considered explore regions of phase space with 
essentially regular classical dynamics when K=0.5, while 
for K=1.5 these regions correspond to a mixed classical 
dynamics. We have also obtained similar results when 
we considered the case with K = 1. 

The range of values of x for which the approximation in 
Ea. (pn|l works is showed in Fig. [71 where we plot the first 
peak of 2?„, that corresponds to the maximum separation 
between the classical and the quantum distributions, as 
a function of x- We can see that the proportionality of 
'Dn with X is in general valid for x < 1 with a slightly 
dependence with the position in phase space of the initial 
distribution. We also observe that, when the phase space 
regions explored by the initial distribution are regular or 
mixed the proportionality is valid up to values I?„ < 1 
while if the regions explored are essentially chaotic the 
proportionality with x is valid up to I?„ w 1 (see inset of 
Fig.EJ. 

It is interesting to compare the behavior of the time 
position of the first peak of 2?^ as a function of ri of our 
curves with its behavior when K = 2 found in |l5j. In 
that case the initial distribution was centered at the ori- 
gin of phase space, that for if = 2 it is an hyperbolic 
fixed point surrounded by a chaotic region (Fig. QJd)). 
The time position of the peak occurred approximately 
for ripeak('7) ~ (l/-'^) ln(l/?7) where A is the expansion 
eigenvalue of the linear map at the origin. However, in 
our case the time position of the first peak of is almost 
the same for all values of rj. Note that the first peak of 
T^n points the moment when the decoherence effects start 
to dominate, restraining the tendency to the separation 
of the classical and quantum distributions produced by 
the presence of nonlinearities in the system. When the 
underlying classical dynamics is chaotic these nonlinear- 
ities are reached in a very short time, so the diffusion 
dominates when these nonlinearities have already been 
reached. This reflects in the behavior of the position of 
the peak that it is a memory of the behavior of the sep- 
aration time Ea. H10|l without the action of a reservoir, 
although lost its meaning due to the smallness of I?„. 
On the contrary, when the underlying classical dynamics 
is regular or mixed the time to reach the nonlinearities is 
very slow so almost immediately the diffusion washed out 
any memory, in the position of peak, of the separation 
time in Ea. ((TTl) . 

Our numerical simulations together with those in [l5l | 
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show that the scale-law in Ea. (|20(l rules the quantum to 
classical transition of the KHO in the semiclassical regime 
and in the presence of decoherence, independently of the 
underlying classical behavior. This is a consequence of 
the approximation in Ea. (|18|l . Therefore, the same scale- 
law proportional to x is valid for the mean values of any 
smooth observables throughout the evolution of the sys- 
tem. Indeed, the mean value of an observable O in a 
state represented by the Wigner function W{n) is given 
by 

(6) = j dxO(x) VK(x) , (21) 

where 0(x) is the Weyl symbol of O obtained through 
the Weyl transform in Eq.Q) replacing p/Amf — > O 3). 
The classical mean value is obtained in the same way 
substituting 0(x) by the classical function O'^'(x) and 
W{x) by the classical distribution P^^'(x). So, for every 
smooth operator in the semiclassical regime we can take 
0(x) « C'(x) |23| and in the presence of decoherence 
we use the result in Ea. lfTB|) to get 




(22) 

V. FINAL REMARKS 

We study the effect of decoherence in the quantum- 
classical transition of the KHO when the classical dy- 
namics is regular or mixed (with small regions of chaotic 
behavior). Our approach is based on the behavior of the 
distance Vn [Ea. H12|l ] as a function of time. We show that 
without the presence of a reservoir, 2)„ clearly indicates 
that the quantum unitary evolution separates from the 
classical evolution because it reaches values that do not 
decrease in the classical limit. When the initial distri- 
bution explore regions of regular classical dynamics, the 
separation time ts scales as a potential law with the ef- 
fective Planck constant, that is, tg « T/h^g. We observe 



that when the degree of mixing increase this potential- 
law is not valid anymore. 

When the system is in contact with a purely diffusive 
reservoir, we extend the results of Ref. 15] for regu- 
lar and mixed classical dynamics, showing that in these 
regimes I?„ is also proportional to x = Kh^f^/^.D^/'^ when 
X < 1. This implies that, in the semiclassical regime 
(?ioff ^1), and independently of the values of the diffu- 
sion constant D and the classical dynamics explored by 
the initial distribution, the quantum-classical differences 
I?„ goes to zero as ^^g- The conditions for the propor- 
tionality between and x can be condensed in the in- 
equality D ^ (if fi^jj)^/'^ ^ hj^g that is the relation of the 
strength of noise and the semiclassical parameter for the 
KHO in order to obtain correctly the classical limit for 
all the regimes of its underlying classical dynamics. Fi- 
nally, we can say that because of the smallness of in 
the semiclassical regime, the concept of separation time 
is not meaningful anymore. 

Although in the framework of an specific model, our 
results and those in R ef. flSl . indicates that the conjec- 
ture presented in Ref. that in the presence of noise 
a single parameter x controls the quantum-classical tran- 
sition of classically chaotic systems, may be extended in 
general for all type of classical dynamics. However, it 
seems that in x should enters only a parameter that con- 
trols the degree of mixing of the considered system and 
not directly the Lyapunov coefficient that characterize 
the exponential divergence of trajectories in fully chaotic 
classical systems. 
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